avances 
en ciencias e 
ingenierias 


Licencia Creative Commons 
Atribucién-NoComercial 4.0 


OOS 


Editado por / 
Edited by: 
Eva O.L. Lantsoght 


Recibido / 
Received: 
06/10/2022 


Aceptado / 
Accepted: 
23 


16/12 


Publicado en linea / 
Published online: 
23/01/2024 


Articulo/Article 
Seccion C/Section C 


Vol. 16, nro.1 
2823 


Numeric simulation and experimental validation for a 
Novel Trapezoidal Solar Collector 


Aien Salvo®™, Pablo Dellicompagni®', Judith Franco®" & Nilsa Sarmiento 
‘Instituto en investigaciones en energias no convencionales (INENCO-CONICET) 
*Consejo nacional en ciencia y técnica (CONICET) 


*Corresponding author: aienweni@gmail.com 


Simulacion numérica y validacion experimental para un 
Colector Solar Trapezoidal novedoso 


Abstract 


In a global warming context, the diversification of the energy matrix is essential for 
mitigation. Solar devices have begun to play an important role in this sense, flat-plate 
solar collectors being the most practical device. On the other hand, numerical modeling 
and experimental validation are important tools for improving the performance of 
these technologies. In this work, a trapezoidal solar air heating collector for food drying 
processes was modeled by using the Simusol open-access software, and experimental 
validation was performed. This particular shape presents a geometrical novelty since no 
other similar designs were found in the available literature, even in such an application 
as food drying. Key parameters, such as air temperature, global efficiency, air mass flow, 
global heat loss coefficient, and useful heat, are determined and discussed. The proposed 
air collector numerically behaves as expected. The output air temperature reaches 
about 100 °C, while the peak heat gain is about 900 W, which makes the air heating 
collector suitable for drying applications. Due to natural convection being the main 
heat transferring mechanism, low air mass flows were obtained. For the case analyzed 
here, this last parameter ranges from 0.012-0.016 kg/s for the optimal thermal behavior. 
The numerical model developed is a reliable tool for designing thermal technologies 
without extra capital cost; the increase of collecting area leads to a considerable increase 
in the thermal power output and improves the psychrometric conditions of heated air. 


Keywords: Solar air heating, trapezoidal collector, computational simulation, experimental 
validation, Simusol software, thermal behavior 


Resumen 

En elactual contexto de calentamiento global, la diversificacidn de la matriz energética es 
fundamental para la mitigacidn. Las tecnologias solares han comenzado a desempenar 
un papel importante en este sentido, siendo los colectores solares de placa plana los 
mas practicos. Por otro lado, el modelado numérico y la validacién experimental son 
herramientas importantes para mejorar el rendimiento de estas tecnologias. En este 
trabajo, se modeld y validé experimentalmente un colector de aire solar trapezoidal para 
procesos de secado de alimentos utilizando el software de cddigo abierto Simusol. Esta 
forma particular se presenta como una novedad geométrica, ya que no se encontraron 
otros disefos similares en la literatura disponible, ni siquiera en aplicaciones como el 
secado de alimentos. Se determinan y discuten pardmetros clave, como la temperatura 
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del aire de salida, la eficiencia global, el flujo de masa de aire, el coeficiente global 
de pérdida de calor y el calor util. Los resultados numéricos del comportamiento del 
colector fueron los esperados. La temperatura del aire de salida alcanza unos 100 °C, 
mientras que la ganancia maxima de calor es de unos 900 W, lo que hace que el colector 
solar sea adecuado para aplicaciones de secado. Debido a que la conveccién natural es 
el principal mecanismo de transferencia de calor, se obtuvieron bajos flujos de masa de 
aire. Para el caso analizado aqui, este Ultimo parametro oscila entre 0,012 y 0,016 kg/s 
para las condiciones dptimas de funcionamiento. El modelo numérico aqui presentado 
resulta en una herramienta fiable para el disefho de tecnologias térmicas sin coste 
adicional de capital. El aumento de la superficie de captacidn conduce a un aumento 
considerable de la potencia térmica de salida, asf como mejoras en las propiedades 
psicrométricas del aire calentado. 


Palabras clave: calentamiento solar de aire, colector trapezoidal, simulacidn computacional, 
validacion experimental, programa Simusol, comportamiento térmico 


BACKGROUND 


Solar energy has been considered the main source of renewable energy in the last 
years [1]. In this sense, technologies and solar devices have begun to play an important 
role in the global warming and climate change context [2]. The implementation of 
solar energy-based technologies, at all scales, promotes a decrease in both fossil fuel 
consumption and greenhouse gas emissions [3,4]. In the field of solar thermal energy, 
flat plate collectors are one of the most widely used technologies due to their easy 
construction and operation [5]. They are heat exchangers that transform solar radiation 
into thermal energy [6,7]. Solar air heating collectors (SAHC) are the most widely used 
devices due to their simple design, although their efficiency is lower than that of solar 
water heating collectors [8]. However, new designs are continuously being proposed to 
increase their thermal performance and obtain more compact and lighter devices [9- 
12]. Since the SAHC are used in multiple applications, they can be assembled in various 
ways in terms of their configuration and constituent materials [13]. 


Numerical modeling combined with experimental validation has been widely 
implemented as a strategy for design optimization of air heating collectors [14-19]. 
However, several authors have studied the SAHC by focusing on theoretical/numerical 
development. Zulkifle [8] evaluated two types of covers—glass and Fresnel lens—on a 
v-groove SAHC and studied the influence of each cover in terms of heat transfer, pressure 
drop, output temperature, and thermal efficiency. The theoretical and numerical 
approaches were established using Excel software. These authors reported higher 
efficiency for a collector covered with Fresnel lens than for those covered with glass, 
71.18 % vs. 54.10 %, respectively. Moreover, the mass flow rate contributed to efficiency 
enhancement. Bhattacharyya [20] determined by numerical simulation the optimum 
number and dimension of rectangular fins in a single-pass SAHC for two different fin 
thicknesses. The design and optimization of the SAHCs were carried out using the 
Computational Fluid Dynamic Model. Qui Lin [21] made a numerical simulation of SAHC 
with a perforated metallic plate by using a ray-tracing solar load model and obtained the 
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temperature and velocity distribution on the plate collector. Hernandez [22] developed 
an analytical model to assess the thermal performance of three different collectors: 
single-pass, double-pass, and double-flow. A robust formulation was elaborated, based 
on the local energy balance, which made it possible to obtain the efficiency plots for 
collectors, as well as the input/output temperatures. Other important parameters, 
such as heat removal factor, convective coefficients, useful heat gain, and global loss 
coefficients, were formulated. The best thermal behavior was revealed for a double-pass 
collector, with an optical efficiency of up to 60 % and a peak output temperature of 46.4 
°C. This behavior is related to better heat transmission when both upper and bottom air 
fluxes are opposite. 


ore advanced numerical simulations can be found in the literature for solar collector 
optimizations. For example, Zhang et al. [23] have developed a neuronal network to 
optimize a water heating solar collector. In this study, nanoparticles were added to 
he heat transfer fluid with different molar fractions and room temperatures, and 
he performance analysis revealed that thermal conductivity and viscosity could be 
predicted with an acceptable accuracy, but neuronal networks need to be trained. 
anoparticles improve the thermal efficiency while the molar fraction increases. The 
numerical optimization is also possible not only for a single solar device. Hajabdollahi 
et al. [24] studied different combinations for water heating solar collector piping. They 
implemented the fast and elitist non-dominated genetic algorithm to find the values 
of optimum design parameters such as bottom isolation, edge isolation, length of 
collector, width of collector, number of tubes, tube diameter, the number of collectors, 
and arrangement of collectors in the network. Fuel cost was also included in the analysis. 
This research concludes that an increase of 10% in the fuel cost yields a decrease in cost 
ratio by 4.75%. Secondly, an increase of 10% in the total solar energy incident on the 
collectors means a cost ratio decrease by 9.69% while efficiency increases by 0.0178%. 
Another research work of solar collectors includes phase change materials (PCM) in their 
optimization analysis. In this line, Abu-Hamdeh et al. [25] compared the energy efficiency 
between nanofluids containing paraffin wax/graphene and paraffin wax/graphene oxide 
phase change materials in flat plate solar collectors. Their results revealed similarities 
with Hajabdollahi et al. [24], in the sense that including nanoparticles, combined with 
PCM, increases thermal conductivity, and thus better performance is observed, with 
even higher values than pure paraffin. 


ost authors validated their theoretical developments under laboratory conditions and 
experimental tests. Reichl [26] performed three-dimensional fluid dynamic simulations for 
SAHC, where the heat mechanisms of conduction, convection, and radiation were modeled 
in great detail. Their model was solved by the ANSYS FLUENT software. The experimental 
setup was carried out at the Test Lab Solar Thermal Systems of the Fraunhofer Institute for 
Solar Energy Systems (ISE) in Freiburg. The modeled collector was validated by comparison 
of the output temperature obtained under laboratory conditions. A high correlation was 
observed, and therefore, the simulation allowed the authors to predict other parameters, 
such as absorber temperature, heat transported to air, efficiency, and thermal losses. 
Badache [27] compared the performance of two types of unglazed thermal collectors 
with computational fluid dynamics through the FLUENT software. In this case, turbulences 
within the collector were considered with the Reynolds—Normalized Group (RNG) k-e 
turbulence model [28-29]. The experimental tests were performed under a laboratory- 
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controlled environment, where the collector had a vertical position and the unglazed 
covering was, in fact, a metallic absorber on which the radiation from lamps impinged. This 
absorber had a certain number of slots, and the air circulation was forced through them. 
Both numerical and experimental results were focused on thermal efficiency, combining 
three mass fluxes with two radiation values and two different plenums. A higher correlation 
of results with a higher mass flow was observed, while radiation and plenum thickness 
variations had no substantial influence on efficiency. 


The validation of numerical models by experimental investigations under real 
solar conditions is more valuable than that conducted in the laboratory [3,30-34]. 
Esmalie [14] developed a complete energy model of a SAHC to predict the air output 
temperature and velocity profiles and then compared the results with measured 
data. This comparison revealed an average difference of 3.5% and 1.5% between 
predicted and measured output air temperature in natural and forced convection, 
respectively. These values were lower than the results obtained in previous studies, 
which included error percentages in a range from 4% to 10%. Mustafa and Mustafa 
[28] performed the numerical simulation of a quarter-circle-shaped SAHC using 
ANSYS-FLUENT, and a mesh model was created with GAMBIT. They considered 
viscous properties and turbulence and implemented the RNG k-e model to predict 
and visualize heat transfer and flow field throughout the SAHC. The validation of 
numerical results was performed with the experimentally measured data, with a 
reasonable correlation when the air flow temperature difference was plotted against 
the collector radius. Hernandez [30] implemented the Gauss-Seidel iterative method 
for equation solving. This method allowed them to design and evaluate the thermo- 
energetic behavior of SAHC with a flat or V-corrugated absorber plate and a porous 
matrix. Higher air temperature and efficiency are possible due to this porous matrix 
incorporation. For the validation of their numerical model, both the inlet and outlet 
air temperatures were measured with K-type thermocouples, and the useful heat gain 
was compared. This comparison demonstrated a high correlation between numerical 
and experimental results for four winter days with maximum differences below 3 °C, 
while the useful heat gain ranged, peaking at 1200 W. 


INTRODUCTION 


The work presented here is supported by an academic project for developing a solar 
dehydrator of tropical fruits and vegetables. The aim of this article is to test how the 
proposed SAHC operates under different sky conditions, by both experimental and 
numerical analyses of the thermal performance, following the described methodology 
in Fig. 1. The solar energy gain takes place in a trapezoidal-shaped SAHC. This particular 
shape presents a geometrical novelty since no other similar designs were found in the 
available literature. In comparison to a rectangular-shaped collector, the air velocity 
increases along the inner channels from the bottom to the top of the collector. On 
he other hand, the construction of a trapezoidal SAHC could require more complex 
echniques than solar collectors available in the market. Furthermore, as the collector 
studied here has bigger dimensions than common ones, its transportability to the 
emplacement place is more difficult. These technical issues make the trapezoidal SAHC 
ess attractive to be proposed, built, and mounted. This particular shape was selected 
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because it was designed to be implemented in tropical regions where environmental 
humidity could be greater than other locations. The additional collecting area also 
improves the psychrometric properties of the air, mainly its moisture absorption capacity 
and specific enthalpy. 


iesleeies 


| and model adjustment Comparison of results 


Figure 1. Scheme of followed methodology for both numerical simulation and experimental validation. 


The SAHC thermal model proposed here is based on natural convection as the main 
heat transfer mechanism for the air heating process. It is developed using a relatively 
new free software, called Simusol [35], to determine its thermal behavior, taking into 
account its geometrical shape, size, and construction aspects. Key parameters, such as 
air temperature, global efficiency, air mass flow, global heat loss coefficient, and useful 
heat, are determined and discussed. The validated model will allow the study of the 
thermal performance of the collector under different climatic conditions. Even more, a 
simulation of a rectangular-shaped SAHC with the same thermal model is presented in 
this work to compare differences in results between both technologies. 


CASE OF STUDY DESCRIPTION 


Fig. 2 shows the analyzed trapezoidal SAHC. The dimension of the collection area is 2.44 
m x 1m wide (bottom length and top sides, respectively) and 1.9 m length. See Fig. 2a. 
The total trapezoidal area of collection is 3.3 m’. The total inner height is 0.1 m, and the 
cross section decreases along the collector length. Therefore, the inlet area of the flow 
is 0.244 m? and the outlet area is 0.1 m?. This configuration influences air circulation by 
increasing air velocity as the area decreases. 


The covering surface is made of 0.006 m thick alveolar polycarbonate, which was 
considered as a homogeneous mass into the physical model, and the internal solar 
reflection effect was dismissed. The external case is made of 0.001 m galvanized steel 
sheet, and its mass was also dismissed in the model. Additionally, the lateral walls and 
the bottom of the collector had a 0.05 m fiberglass layer. See Fig. 2b. 
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Figure 2. (A) Scheme of the trapezoidal SAHC and its external dimensions in meters. (B) Inside view of the solar 
collector output, to be connected into a drying chamber. 


The absorber plate is 1 mm-thick, black-painted sinusoidal galvanized steel. However, the 
absorber was considered as a plane plate, and its thermal conductivity was neglected in 
the model. The disposition of the absorber plate leads the air circulation into upper and 
bottom channels. The assumed thermal properties of the materials used in this work are 
shown in Table 1. 


Table 1. Properties of the materials 


Property Absorber Covering Case Isolation 
Density ep (kg/m?) 7.850 1.200 > 14 
Conductivity \ (W/m °C) - 0.2 - 0.03 
Absorptivity a (-) 0.96 = = = 
Transmittance t (-) : 0.8 x i 
Emissivity € (-) 08 ty) 0.1 = 


PHYSICAL AND MATHEMATICAL MODEL 


The following hypotheses were considered to formulate the model: 


+ The air circulation inside the channels is unidirectional, in the input-output 
direction. In practice, it is favored by the sinusoidal shape of the absorber. 


+ The roughness of the internal walls is neglected. Under this hypothesis, a uniform 
air velocity profile for each cross section was assumed. 


+ — Transient disturbances of external wind velocity are not considered. 


+ The internal multiple reflection effect inside the covering plate is neglected here 
since there are no precise models related with honeycomb polycarbonates to 
represent this effect in the literature. 
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The thermo-physic properties (thermal conductance, absorbance, and reflectance) 
are considered as unchanging and the composition of materials homogeneous. 
All the heat fluxes are considered unidirectionally. 


The psychrometric properties of the air nodes (absolute and relative humidity) are 
not considered at this stage of simulation. 


The spectral solar transmittance of the polycarbonate covering layer is considered 
as a constant value of 0.85, and its dependence on solar incidence angle is 
dismissed. 


thermal energy is transferred from the collecting plate to the circulating air by heat 
vection as the main mechanism [6]. Then, radiative and conductive transmissions 
also considered. The air circulation takes place both below and above the 


absorber plate. This configuration facilitates improvements in the performance of 


the 


heat collecting process [22,36]. Hourly measured climatological variables, such 


as wind speed, solar radiation, and environmental temperature, were considered as 
boundary conditions. The outer convective coefficients were calculated based on the 
methodology described by Duffie and Beckman [6]. Fig. 3 shows the physical model 
implemented in the simulation. In this model, the air temperatures are represented 
by nodes: one of them located between the plate and the covering surface, and the 
other one between the plate and the bottom surface. The transmitted solar radiation 


(It) 


impinges on the absorber, and part of it is absorbed as thermal energy (/ta), 


which is partially transferred to the air nodes by convection. The thermal transmission 
between these air nodes and lateral inner walls is also considered. Thermal conductive 
mechanisms through the isolation case and the cover are taken into account. Three 
thermal losses take place in this model: upper (covering surface), lower (bottom 


surface), and lateral. The boundary conditions on the outer and inner surfaces also 


include the effect of thermal radiation and convection. 


h2=0.033 m 


Figure 3. Scheme of temperature nodes and thermal transmission circuits, for an arbitrary transversal section 
considered, where dimensions are in meters. The colors for the different heat transmission mechanisms are: red 
for irradiative, blue for convective, and green for conductive. 


DOI: https://doi.org/10.18272/aci.v16i1.2823 


Articulo/Article 
Seccin C/Section C 


Vol. 16, nro. 1 
2823 


avances 
en ciencias e 
Ingenierias 


8 


Numeric simulation and experimental validation for a Novel Trapezoidal Solar Collector 
Salvo / Dellicompagni / Franco / Sarmiento (2024) 


The ribbed shape of the plate is not considered in the model because of its quasi-flat 
shape. A different case occurs when the absorber has a pronounced irregular form, 
as noted by Zulkifle [8] and Bashria [37]. Nusselt number conditions for laminar flow, 
transition, and turbulent flow regions must be treated differently. Nevertheless, in the 
present design, the transversal areas where air flows can be considered rectangular, 
without much error. The equations for thermal calculus are described below [38,39]. 


_ pi Vi Dhij (1) 
i hi 
_ 0.453 Re;;°° : Pr Pri; > 0.6 
4 (0.0308 Re; ;°% - Pr;;°393 (0.6 < Prij < 60 2) 
2 gg 3 
hy = Shy (3) 


The temperature dependence of Prandtl (Pr), dynamic viscosity (u), and conductivity 
coefficient (A) are implemented within the program by tables. The values for these 
parameters were taken from Bergman (2011). A similar consideration was reported by 
Zulkifle [8] and Bashria [37], but they implemented formulas related to air temperature, 
instead of discrete values. 


As the collector was modeled by dividing the equipment into three sections, subscript / 
refers to each transversal discrete section of the collector while subscript / refers to each 
convective and radiative coefficient within the same section. Each section is connected 
to the others by the air flow mass, as shown in Fig. 4. B/-3 and D1-3 nodes represent the 
temperature of circulating air above and below the absorber plate, respectively. The air 
mass flow through the collector is the same, but since the transverse section has a linear 
decrease, the air velocity increases while the air rises through the collector. H1-3 and L 1-3 
represent the inertial capacitance of the circulating air. The temperature nodes ss and si 
represent the output air temperature for both upper and bottom channels, respectively. 
Convection was modeled between these nodes and the environmental temperature 
node, according to Iriarte [40]. 


DOI: https://doi.org/10.18272/aci.v16i1.2823 


Articulo/Article 


Secci6n C/Section C 
Vol. 16, nro. 1 Numeric simulation and experimental validation for a Novel Trapezoidal Solar Collector 
2823 Salvo / Dellicompagni / Franco / Sarmiento (2024) 
avances Air flow Air temperature node 
encciencias e 


ingenierias 


Figure 4. Model representation of air mass flow through the collector. BO-3 and DO-3 correspond to the air flow 
above and below the absorber plate, respectively. 


Boundary conditions 


For any section, the mathematical boundary conditions for the calculation domain are 
given as: 


* Covering surface 
This surface is exposed to the outer conditions. Convection and radiative 
transmissions are the heat transfer mechanisms. The aperture area is made of 
honeycomb polycarbonate and modeled as single thermal conductance with a 
transmission coefficient of A,=3.6 (W/m? °C). Eq. 4 gives the energy balance for 
the solar collecting surface. 


T ii-T ci 
Ns,i (Tci,o a Te) + FstyBs€cO(Tei,o* = Tsky*) te Fspy qa ~ Bs)éco(Tei,o* — Te*) = Ap fetal (4) 


According to Liu [40-43]: 


B, = §0.5(1 + cos 6;) (5) 
Fsxy = 0.5(1 + cos 6;) (6) 
Terry = 0.05527,” (7) 


- Absorber plate 
The absorber is a black painted metallic surface, where solar radiation is partially 
absorbed and the temperature of the plate rises. The energy balance is given by 
Eq. 8, where thermal conduction through the metallic surface was neglected due 
to its thinness. 


[Tp Qp _ Apai,a (1, \ a Taia) + Apaip(T; i Tai,b) a €,0[2Tpi* = Tit? _ Tei’ | (8) 
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- Bottom surface 
The air circulation below the absorber transmits part of its thermal energy to the 
bottom surface by convection, and then by conduction through the bottom 
isolation layer to the outer surface. For any section of the collector, Eq. 9 introduces 
the thermal condition boundary for the bottom surface, where the conduction of 
the metallic internal case is neglected due to its thinness. 


T pii-T bio 
hav i(Taib — Toit) + €0(Tpi* — Trii*) = Ais Coss tors) ©) 


- Lateral inner and outer walls 
The thermal losses through the lateral walls of the collector were also considered 
in this model. Only the convective mechanism is considered here, for both bottom 
and upper channels: 


h [ees (ig = Tae 4a li,i-2 Gag =F eds (Tii,i-1-T to i-1) (10) 


ai,a ai,b Ax, 
lii-2 lii-2 _ 4. (ui-2-Ttoi-2) (11) 
AV ig (Tota — Tug-2) +h |r, (Tote —Titi-2) = As he 


Then, the heat throughout the walls is transferred to the external atmosphere by 
convective and radiative mechanisms (Eqs. 12 and 13). 


qe Vieni) a | was Con _ Te) + d(T" = 7?) (12) 


Nis Pesce tate) = eee Cen = Te) + €0(Tio,i-2" = Te) 


Thermal performance 


By definition [6], the thermal loss represents the energy that is not transferred to the 
heat carrier fluid (HTF), and this energy goes to the atmosphere. A global thermal loss 
coefficient is defined in terms of convective, conductive, and radiative coefficients. It is 
also possible to calculate the thermal losses through the top, bottom, and lateral walls 
separately (Eqs. 14, 15, and 16, respectively). 


Ur = Di} Oa,,i [Aci (Tei,i — i)\ (14) 

Us, S03 9,,.,i|Abi (Toi = Taal} (15) 
j=2 

U, = DIF Oy,,1-j|Aui—-j(Tui-j — Tioi-;)] , (16) 
j=l 


Heat fluxes Dr, i D,,,, ANG Dy,.i; are calculated by Simusol as conductive losses, and their 
values can be obtained by simply requesting the simulation results from the software 
program. Eq. 17 gives the global heat loss coefficient [6]. 


Ug =U, + Up, + U, (17) 
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The collector's useful energy gain and the global efficiency are calculated by means of 
Eqs. 18 and 19, respectively. 


Qu = Cp[tta (Taiali-g — Te) + Mb(Taipl-3 — Te)] (18) 


Ng = 100 « # (19) 


Cc 


Both air-mass flows m, and m, are determined by instantaneous thermal balance at air 
temperature nodes B/ and D1], respectively [40]. 


j=4 


Mg Cy Taialiey = j=l Da jlint (20) 
My Cy Toi plo. = ye Dp flint (21) 

So, the total air mass flow results in: 
Tit, = Tg + Thy (22) 


The convective fluxes Dajfie1 aNd Byjfi-1 are calculated by Simusol and represent the 
thermal gain from the absorber, lateral walls, covering, and bottom surface into 
respective air flow nodes, as Fig. 5 shows. 


inner cover surface 
upper 


an convective flux 


absorber 


inner lateral surface inner lateral surface 


hatin convective flux 


air . 
inner bottom surface 


Figure 5. Scheme of convective fluxes into air flow nodes, for the first section of the SAHC 


SIMULATION RESULTS 


Simusol generates a spreadsheet with a 7 s time-step for each parameter. It also 
shows a set of time-dependent variables upon request; otherwise, Simusol shows the 
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temperature profile for all the thermal nodes in the model. To perform the simulation, 
the model presented here requires some climatological data input: the solar radiation 
based on experimental measurements, the environmental temperature, and the wind 
velocity in the area surrounding the collector. This information is given by tables, and 
they can be repeated many times by a Prep command. It is also important to define the 
slope of the collector and the geometrical dimensions. Since Simusol runs on Linux, it 
can be launched in a user terminal from the folder where the .dia file is located, using the 
simusol command. It is important to save all changes made in the model before starting 
the simulation. Fig. 6 shows a diagram of the interaction between the thermal model, 
the climatological data entry, the thermal/optical properties of the materials, and the 
geometrical design parameters for upper and bottom channels. Differential equations 
are solved by the fourth-order Runge-Kutta method [41], which is a simple, robust, and 
suitable option for the numerical solution of differential equations combined with an 
intelligent adaptive step-size routine. 


INPUT PARAMETERS 


Climatic variables: 
~ Solar radiation 
~ Environment temperature 


OUTPUT PARAMETERS 
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Figure 6. Scheme of the thermal model implemented on Simusol. 


Four days were selected to perform the numerical simulation, with mainly two types of 
days: clear and partially clouded. The aim was to reproduce the conditions of variability, 
as realistically as possible, with which the experimental SAHC must operate. Fig. 7 shows 
the temperature profile of the air flux in the upper channel. To simplify the graphs, a 
period from 9 a.m. to 5 p.m. was defined as the operating time of the SAHC. 
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Figure 7. Simulated air temperature for the upper channel. Temp 1 to 3 and Temp out correspond to nodes B 1 
to 3 and temp ss in Fig. 3. (A) March 29, (B) April 6, (C) April 12, (D) April 15. 


It can be observed that the output temperature of the upper air can reach about 100 °C. 
For some processes like food drying, this temperature could be higher than necessary but 
ideal for certain grains, seeds, and corn. However, this temperature value enables the SAHC 
to be applied to other processes with higher thermal demand. The high temperatures 
obtained here are strongly related to the high solar radiation values measured at the 
location and the additional collecting area due to the trapezoidal shape of the SAHC. 
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Figure 8. Air mass flow profiles for both upper and bottom channels of the SAHC and total air mass flow. 
(A) March 29, (B) April 6, (C) April 12, (D) April 15. 
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Eqs. 20, 21, and 22 give the air mass flow in kg/s for both upper and bottom channels and 
the total air mass flow, respectively. The heat gain and the air temperature are important 
in establishing the air mass flow (Fig. 8). It can be observed that the air mass flow is 
different for each channel because of the different heat gains. The highest air mass flow 
corresponds to the upper channel, which means that the air velocity is the highest there. 
This numerical result was expected because the heat source (the absorber) is located 
under the air thermal node in the upper channel, while the same heat source is located 
above the air thermal node in the bottom channel. 


The values of the total air mass flow range from 0.0015 kg/s to 0.02 kg/s. Other studies of 
SAHC yielded values of 0.019-0.028 kg/s [20], 0.017 kg/s [44], 0.013-0.04 kg/s [45], and 
0.056 kg/s [46]. It is important to note that there are significant differences between the 
mentioned systems, such as the position and slope, the geometry, and the absorber 
shape and materials. Nevertheless, the air mass flows obtained here correspond to those 
normally obtained for this kind of technology. 


The useful heat gain is determined by Eq. 18, which takes into account the air mass flows of 
both channels and the difference in output-input air temperatures. Here, the dependence 
on the temperature of the specific heat of the air was neglected. Fig. 9 shows the useful 
heat profiles for the four days considered. A maximum thermal power output of about 
1200 W can be reached by the SAHC under study. The global efficiency is determined 
by Eq. 19, and Fig. 9 shows this parameter, where, in clear-sky conditions of operation, an 
average value of 35% is observed, in the hours around the midday, while efficiency varies 
on cloudy-sky days. Peaks of about 45% on global efficiency are reached by this SAHC. 
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Figure 9. (A) Global efficiency, and (B) useful heat gain. 


One of the main parameters to be considered in describing the collector's performance 
is the global heat loss coefficient, which is given by Eq. 17. In this model, the top, lateral, 
and bottom thermal losses were considered separately, and then the global heat loss 
was estimated. Fig. 10 shows the time variation of the total global heat loss coefficients 
for all the considered days. For the operating period, the global heat coefficient ranged 
from 16 to 18 W/m?*°C. The literature reports a normal value of 10 W/m?’C [6] for a single 
covered air heating solar collector. 
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Figure 10. (A) Global heat loss coefficients. (B) Useful heat gain related to the total air mass flow. 


Fig. 10 shows the range of air mass flow for the heated air production, considering the 
thermal power output of the SAHC. For optimal performance, the SAHC should operate 
in a range of 0.016-0.020 kg/s. This must be considered if a forced air circulation system 
is integrated into the SAHC. 


EXPERIMENTAL SETUP 


The trapezoidal SAHC location was -24.69 Latitude and -65.39 Longitude, and it was 
north-oriented with a fixed optimal slope of 36° for a homogeneous solar collection 
throughout the year. In order to replicate the real working conditions, the SAHC was 
connected to a drying chamber of 1 m? volume, as shown in Fig. 11a. 


Fluke 
Thermometer 


Hot Wire 
_ Anemometer 


Solar 
| Meter 


Figure 11. (A) Trapezoidal SAHC connected to the drying chamber, in its working position. (B) Scheme of 
measuring devices used for data recording. 


The solar radiation was measured with a solarimeter placed on the covering surface of the 
SAHC, as shown in Fig. 11.b. The air temperatures were measured according to the nodes 
in the simulation model: B1, B2, B3, and ss. It was only possible to place the temperature 
sensors in the upper channel of the collector. A thermometer Fluke 54 Il with type-K 
thermocouples was used for measuring the temperature. The surrounding wind velocity 
was measured using a hot-wire anemometer. Table 2 briefly describes the used instruments. 
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Table 2. Properties of the measurement equipment 


Parameter Device Accuracy 

Temperature Fluke 54 Il thermometer with thermocouple K-type 0.1 °C 

Wind velocity Hot-wire anemometer TES 1341 0.01 m/s 

Solar radiation HT204 solar power meter 1 W/m2 
Experimental results 


The temperatures of the circulating upper air were measured. As commented in the 
previous section, K-type sensors were used and the measuring procedure was carried 
out during the operative period of the SAHC, from 9 a.m. to 5 p.m. Although the ASHRAE 
code sets the conditions for testing the equipment [47], both clear-sky and cloudy-sky 
conditions were considered to make simulations under real operative situations. Fig. 12 
shows the measured air temperatures for the upper channel. It can be observed that 
the air temperature peaks between 80 °C and 100 °C. As expected, the shape of the 
curves follows the solar radiation profile. This indicates that the SAHC has low thermal 
inertia, and it can quickly reach the thermal state of operation when the solar radiation 
has suitable values. Since drying processes require a temperature range of 50-60 °C [48] 
for optimal performance, the SAHC needs solar radiation between 600 W/m? and 900 
W/m?. For cloudy-sky conditions, with solar radiation of about 200 W/m”, the outer air 
temperature is about 35-40 °C. This represents an advantage because a reduced size of 
an auxiliary heat source is available, if necessary. 
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Figure 12. Measured air temperatures in the upper channel. The Jemp out plot corresponds to the ss node from 
the simulation. (A) March 29, (B) April 6, (C) April 12, (D) April 15. 
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The main reason for these differences is the weather input data values for every 
single day, referred basically to the solar radiation, environmental temperature, and 
surrounding wind velocity. It can be observed that the temperature response in the 
simulations coincides with the experimental measured one (every single plot in Fig. 7 
compared with those in Fig. 12). 


MODEL VALIDATION 


The experimental and numerical results were compared (Fig. 13). Based on direct 
measurements, the air temperature for the upper channel was considered the main 
variable for performing the comparison. For the set of measured days, the determination 
coefficient R* ranged from 0.82 to 0.93. The lowest correlation was obtained for a mostly 
cloudy day, where the predominant solar radiation was about 200 W/m? after midday 
(April 12). This correlation can be considered acceptable despite the unstable working 
conditions. For a clear-sky condition, the correlation coefficient was 93%. Intermediate 
correlation values were observed for semi-cloudy sky conditions. The more stable 
climatic conditions, the better adjustment is achieved by the numerical model. 
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Figure 13. Comparison between experimental and numerical air temperatures. 
(A) March 29, (B) April 6, (C) April 12, (D) April 15. 


The scattering observed in Fig. 13 can also be explained by the surrounding wind 
variability when the measurements were performed (Fig. 14). As the SAHC studied here 
is a natural convective one, the oscillation of external wind produces variations in the 
external pressure near the inlet/outlet zones, influencing the air mass flow inside the 
SAHC, and therefore, the internal convective coefficients present high variability. Fig. 14 
also shows a range of 0.1-3.4 m/s for wind velocity. This is a climatological feature of the 
region where the SAHC is mounted [49]. 
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Figure 14. External wind velocity surrounding the SAHC. 


The reliability of the comparison was analyzed using Infostat software. In addition to the 
R*, other coefficients were determined. Table 3 shows the p-value, standard error (Gg), 

x 
mean absolute error (€), asymmetry, and Kurtosis coefficients. 


Table 3. Statistical parameters for the model validation process 


Asymmetry Kurtosis 
Meas. Sim. Meas. Sim. 
29% March 92.81 0.0000 5.86 4.80 1.14 0.44 -1.11 -1.78 
6" April 87.36 0.0000 7.2) Soy 1.89 0.68 -0.84 1.57 
12" April 82.36 0.0000 5.78 443 6.63 455 5.95 2.39 
15% April 84.14 0.0000 9.69 8.15 (33) 0.34 = 11.55 = 125) 
The results of the statistical analysis show a p-value lower than 0.05 for all cases, which 


indicates a statistically significant relation between the simulated and measured values, 
with a confidence level of 95%. As mentioned above, the model explains the correlation 
between the analyzed variables by an R? higher than 80%, which indicates a good 
relation between variables. In general, asymmetry and kurtosis coefficients present 
values within a 2 — -2 range. This implies that the data sets have a normal distribution. 
However, it was observed outside the range values only for April 12. This significant 
deviation was caused by the high meteorological variation: solar radiation decrement 
combined with high external wind velocity. 


DISCUSSION 


The numerical results obtained in the present study are closely comparable with the results 
of other authors. The comparison of the main parameters is shown in Table 4. As mentioned 
in the introduction section, no one single case for trapezoidal collector was found in the 
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literature. While most of the designs were focused on thermal efficiency, with values 
above 60% (similar to vacuum pipes based-on collector) mostly reached, low thermal 
power outputs were reported. Additionally, the reported outlet air temperature was no 
higher than 71 °C. This could be related to different parameters such as solar resource and 
sky conditions, orientation and disposition of collectors, sizes of solar absorber, materials 
used for prototype constructions, etc. The solar collector studied here reached higher 
temperatures and power thermal output, despite the fact that lower thermal efficiency 
was performed, compared with other works. This increase of air temperature and power 
thermal output was due to the extra area conforming the trapezoidal shape of the 
collector. Additionally, a numerical scenario as the SAHC built with rectangular shape was 
presented here. In this case, temperatures, efficiencies, and power outputs were similar to 
other authors. Furthermore, the comparison between both proposed models (trapezoidal 
and rectangular) revealed that better air heating conditions are easily achieved by simple 
modification in geometrical shape and collecting area expansion. 


Table 4. Comparison between results obtained in other studies 


, ; Useful 
Author Technology Air mass flow Tout Efficiency eet 
Kareem et | Forced multi-pass Q 
al, 2016 solar dryer 0.01-0.07 kg/s 175m 55.49 °C 59.96 % = 
Lietal Vertical SAHC 
e with aluminum 0.03 kg/s 04m? | 164°Cav. 80-89% | 202Wav. 
2011 
perforated plate 
aaeseel V-shaped flat-plate 
3017 "| SAHC with micro- 0.054 kg/m?s 1.8m? 465°C 69% = 
heat pipe arrays 
ia \V-grove with: ) 35,85¢C 47.10% 
Zulkifle et : ; a) 32: a) 47.10% 
al,201g | ®singleglass ode "| p)3839ec |b) 71.18% : 
b) Fresnel lens 
Present Trapezoidal 0.0015-0.02 ka/s | 3268m? | 1035°C 35% av. 1321W 
study shaped SAHC : id ‘ , 45% max. max. 
Rectangular 21% av. 
aac . 0,0015-0.013 ko/s | 1.9m? 72°C : tg 
study shaped SAHC 49.7% max. max. 
CONCLUSION 


A numerical model is proposed based mainly on a convective heat transfer mechanism. 
The air heat gain was modeled for both upper and bottom channels. This model differs 
from those of other authors in that the air temperatures are considered individual thermal 
nodes, so it is possible to determine each heat flux in the air nodes. Modeling based on 
convection as the predominant heat transfer mechanism made it possible to obtain 
reasonable numerical results and determine the most suitable materials for air heating 
applications. The numerical analysis of the trapezoidal-shaped SAHC considers the 
measured values of the solar irradiation that reaches the outer covering surface, as well as 
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the environmental temperature and external wind velocity. Two possible scenarios were 
considered to test the model proposed: clear-sky and cloudy-sky conditions. The numerical 
simulation was performed using an open-access software program called Simusol. 


It was found that the SAHC can heat the air up to 100 °C. Moreover, lower temperatures can 
be obtained from lower solar radiation, of over 50 °C. The air mass flow ranged between 
0.0015 and 0.02 kg/s, with an optimal range of 0.016-0.02 kg/s for the maximum heat gain 
(1140-1321 W). The numerical simulations revealed a satisfactory thermal performance 
for food drying applications and other low enthalpy processes. It is important to note that 
the trapezoidal-shaped SAHC studied here is a natural convective type and the obtained 
performance has advantages over similar technologies since an air forcing system is not 
needed, allowing energy savings. Although the achieved global efficiencies are not the 
main feature of the proposed design, 35% under stable operating conditions would be 
considered acceptable if the thermal performance (output temperatures, power thermal 
generation, and global heat loss) is taken into account. 


The trapezoidal design of the collector increases the collecting area in a range of 1.8 
to 8.2 times higher than the rectangular technologies found in the literature, without 
extra capital cost or initial investment and using the same constructive techniques. By 
increasing the collecting area, a considerable increase in the thermal power output was 
noted, despite the decrease in the global efficiency. Furthermore, Table 4 shows the 
numerical results for a rectangular SAHC with 1.9 m? of collecting area (1 m width and 
.9 m length), which was simulated by using the validated model presented here and 
taking in consideration both the same climatic conditions and constructive materials as 
the trapezoidal one. It is observed that the output temperature remains about 100 °C 
despite its poor global efficiency. Low useful heat is also observed for this rectangular 
configuration, and this could be explained due to the low air mas flow under natural 
convective conditions, suggesting the implementation of air forced systems to the 
designed SAHC in order to improve the thermal gain and the efficiency as well. 


The modeling implemented in Simusol demonstrates once again that it is a reliable 
numerical tool for designing thermal technologies. Although the SAHC behaved as 
expected, more experiments are necessary to fine-tune the model. Additionally, a more 
accurate numerical description of the internal convective heat transfer mechanism 
could further improve the model. More accurate heat and mass transfer modeling will 
be done in the future, even with the drying chamber included in the model. 


NOMENCLATURE 


Ay; Area of the bottom layer, for the / section, m? 

A. Total area of the covering surface, m? 

Ai;; Area of the lateral layer, for the / section andj element, m? 

A:; | Area of the covering layer, for the i section, m? 

C, Average heat capacity of the air, J/kg °C 

Dhj Hydraulic diameter of channels, for the j section and j element, m 

Fy — View sky factor, dimensionless 

ha; Convective coefficient between the bottom surface and bottom air nodes, for the 
i section, W/m?2°C 


DOI: https://doi.org/10.18272/aci.v16i1.2823 


Articulo/Article 


Secci6n C/Section C 
Vol. 16, nro. 1 Numeric simulation and experimental validation for a Novel Trapezoidal Solar Collector 
2823 Salvo / Dellicompagni / Franco / Sarmiento (2024) 


A Convective coefficient from outer cover surface to the environmental nodes, for 
thei section, W/m?°C 
avanees. Internal convective coefficient, for the i section and j element, W/m?2°C 
mene pli. Convective coefficient between the upper air node and left inner wall, for the i=/ 
section, W/m?°C 
All? Convective coefficient between the upper air node and right inner wall, for the i=7 
section, W/m?°C 
hl!2,, Convective coefficient between the bottom air node and left inner wall, for the j=7 
section, W/m?2°C 
hl’? Convective coefficient between the bottom air node and right inner wall, for the 
i=1 section, W/m?°C 
h Convective coefficient between the absorber and the upper air nodes, for the / 
section, W/m?2°C 
Convective coefficient between the absorber and the bottom air nodes, for the / 
section, W/m?°C 
l,l Solar radiation that impinges into the covering surface, W/m? 
Upper air mass flow, kg/s 
m, Bottom air mass flow, kg/s 
m, — Total air mass flow, kg/s 
Nu usselt number, for the jsection and j element, dimensionless 
Pr, Prandtl number, for the for the / section and j element, dimensionless 
Qu = SAHC useful energy gain of the output air, W 


Re, Reynolds number, for the j section and j element, dimensionless 


Taig Temperature of the upper air node, for the i section, °C 
Tiial)-1 Fe Mperature of the upper air node, for the i=1 section, °C 
Tgicli-3 Temperature of the upper air node, for the /=3 section, °C 

T,» Temperature of the bottom air node, for the i section, °C 
Taiol)-1 Fe Mperature of the bottom air node, for the i=7 section, °C 
Tipl)-3 Temperature of the bottom air node, for the i=3 section, °C 
T,,, Temperature of the inner bottom surface, for the / section, °C 
Ty. Temperature of the outer bottom surface, for the / section, °C 
T. Temperature of covering layer, °C 

T,, | Temperature of inner cover surface, for the i section, °C 

Tj. Temperature of outer cover surface, for the i section, °C 

T, Environmental temperature, °C 

T, Temperature of bottom insulation layer, °C 

Tj; Temperature of left lateral inner wall, for the /=7 section, °C 
Tj2 Temperature of right lateral inner wall, for the i=7 section, °C 
T,,; _ Temperature of lateral inner wall, for the / section and j element, °C 
T,., Temperature of left insulation layer, °C 

T.2 Temperature of right insulation layer, °C 

Ti:, Temperature of lateral outer wall, for the isection and j element, °C 
T,, Temperature of the absorber, for the i section, °C 

Ty, Temperature of sky, °C 

U, Total heat loss coefficient for the bottom layer, W/m?°C 

U, Global heat loss coefficient for SAHC, W/m? °C 

U, — Total heat loss coefficient for the lateral layers, W/m?°C 

U, — Total heat loss coefficient for the covering layer, W/m?°C 
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V, Internal air velocity, m/s 


i 


Y Thickness of the covering/bottom/lateral layers for conductive fluxes calculation, m 


GREEK SYMBOLS 


@,,,; Conductive heat flux throughout the bottom layer, for the / section, W 

@,,,,i; Conductive heat flux throughout the lateral layer, for the i section and j element, W 

Dy, Conductive heat flux throughout the covering layer, for the i section, W 

Dyhi=1 Convective heat flux from jinner surface to the upper air node, for the i=7 section, W 

Dyji-1 Convective heat flux from j inner surface to the bottom air node, for the i=1 
section, W 

a,a Absorptance of the metallic absorber, dimensionless 

B Factor that splits the heat exchange with the sky dome between sky and air 
radiation, dimensionless 

Nz Global efficiency of the SAHC, % 

6, Solar incidence angle, ° 

A Conductive coefficient of materials, W/m °C 

A; — Air conductivity coefficient, for the i section and j element, W/m °C 

Xi Conductive coefficient of fiberglass, W/m °C 

A, Conductive coefficient of polycarbonate, W/m °C 

uw; ~~ Dynamic viscosity, Pa s 

TT Transmittance of the polycarbonate covering layer, dimensionless 

€ Emissivity coefficient, dimensionless 

p Density of materials, kg/m? 

p _ Air density, kg/m? 

oO Stefan-Boltzmann constant, 5.670373 x 10° W/s K* 
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